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Abstract 

We investigate the propagation of a charged particle in a spatially constant, but time de- 
pendent, pseudoscalar background. Physically this pseudoscalar background could be provided 
by a relic axion density. The background leads to an explicit breaking of Lorentz invariance; 
as a consequence the process p — s- p7 is possible and the background acts as a shield against 
extremely energetic cosmic rays, an effect somewhat similar to the GZK cut-off effect. The 
effect is model independent and can be computed exactly. The hypothetical detection of the 
photons radiated via this mechanism would provide an indirect way of verifying the cosmological 
relevance of axions. 
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1 Introduction 



Cold relic axions resulting from vacuum misalignment [U [2] in the early universe is a popular and 
so far viable candidate to dark matter. If we assume that cold axions are the only contributors to 
the matter density of the universe apart from ordinary baryonic matter its density must be [3] 

p ~ 10"3°gcm-3 ~ IQ-^^GeV^ (1) 

Of course dark matter is not uniformly distributed, its distribution traces that of visible matter 
(or rather the other way round). The galactic halo of dark matter (assumed to consist of axions) 
would correspond to a typical value for the density [1] 

Pa ^ lO^^^gcm-^ ~ 10-^°GeV^ (2) 

extending over a distance of 30 to 100 kpc in a galaxy such as the Milky Way. Precise details of 
the density profile are not so important at this point. Even in this case this is a rather low density; 
the axion background provides a very diffuse concentration of pseudoscalar particles, but then the 
photon density is also low and yet it has an impact on ultra-high energy cosmic rays imposing the 
GZK cutoff [5]. What are the consequences of this difuse axion background on high-energy cosmic 
ray propagation? This is the question we answer in this work. 

The fact that the axion is a pseudoscalar, being the pseudo Goldstone boson of the broken 
Peccei-Quinn symmetry [6j, is quite relevant. Its coupling to photons will take place through the 
anomaly term; hence the coefficient is easily calculable once the axion model is known 

A/: = ga^^ — —FF. (3) 

Two popular (and so far viable) axion models are the DFSZ[8j and the KSVZfT] ones . In both 
models ga-y-y — 1. Here a is the axion field while fa, the axion decay constant, is bound to be of 
0(10^^) GeV or more. As to the axion mass, rua it should lie in the range 10"'^ eV > nia > 10~^ 
eV. 

The axion field is originally misaligned and in the process of relaxing to the equilibrium configu- 
ration coherent oscillations with q = are produced, provided that the reheating temperature after 
inflation is below the Peccei-Quinn transition scale. In late times the axion field evolves according 
to 

a{t) = gq cos mat. (4) 
Then, integrating by parts in the previous equation we get 

AC = -ga^^-^sm{mat) e'^^AiFjt (5) 

J a 
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A particle travelling at the speed of light will see coherent regions with quasi-constant values of the 
axion background of a size ranging from the milimeter to the meter, depending on the axion mass. 
Small as this size is, it is very many orders of magnitude bigger that the wave length of a particle 
travelling with a momentum p characteristic of a very high energy cosmic ray. For processes with 
a high momentum transfer we can treat this slowly varying term as a constant and thus replace 
© by 

A£ = ^r^e'^'^A.F.k, (6) 

where the "constant" rj will change sign with a period l/nia- r] is obviously related to the axion 
density. Numerically we expect 

2a ^/p^ ,„_24 in-25. 



M ga-y, — ^ 10-^* - 10-^^eV, (7) 

la 

for fa in the range 10^^ — 10^^ GeV and where we have included and additional factor of ^ to 
roughly account for the variability of the sin mat term. Equation ([6|) can also be written as 

A£ = ^r^cApF'^^ (8) 

with r/o, = 0, 0, 0). The expression uses a covariant language but it is of course Lorentz invariance 
violating as rja is a constant axial vector. 

2 Solving QED in a cold axion background 

Let us study the effects of explicit breaking of Lorentz invariance by means of a time-like constant 
axial vector. Consider electromagnetism in such a background 

= C,my + C,iA\ (9) 

It will be useful for us to keep niy > 0. It need not be a fundamental mass; but an effective or 
induced in-medium mass. Otherwise gauge invariance is manifest. 
The equations of motion are 

{g^-^ {e - + ie^-^"" Va k^} Ax{k) = 0. (11) 

We can build two complex and space-like chiral polarization vectors e^^(/c) which satisfy the 
orthonormality relations 

- g,, e^*{k) el{k) = 1 g e^*{k) e^k) = (12) 
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In addition we have 

£!^(k) ~ F eHk) ~ k'^T]'' -k^'r]-k (13) 

They fulfih 

e^;{k)e-B{k) = gAB ^^(fe) e'^B W = 9^"- (14) 

These polarizations correspond only approximately to the usual ones of QED. As a consequence 
light propagation in an axion background may be subject to modifications. However, for visible 
light and even less for radiowaves, there is no marked separation of scales with respect to the axion 
mass and therefore the time variation of the background cannot be treated adiabatically. The net 
effect is expected to be zero or extremely small [IJJ. We refer the reader to |9i |10j for additional 
details. 

The polarization vectors of positive and negative chirality are solutions of the vector field 
equations if and only if 

A;^ = (cuk±,k) a;k± = ^k2+m2±r/|k| e^(k, t?) = 4(A:±) ( 4 = c^k± ) • (15) 

In order to avoid problems with causality we want k'j_. > 0. For photons of a given chirality (negative 
if ?7 > 0, positive if r/ < 0) this can be if and only if 

2 

TTl 

|k|<^ = A^. (16) 

7] 

In fact for = these photons cannot exist as physical asymptotic states. If they are produced 
they will eventually decay (to three photons of like chirality) in a cascade process that leads to a 
red-shift. 

As is known to everyone the processes e~ e~7 or 7 ^ e~^e~ cannot occur in vacuum. However 
in the present situation, because of the preferred reference frame, physics is different in different 
frames and for the latter process 

^^k± = \J^'^ + m^, lb r/ |k| = \/p2 + + y^(p - k)2 + m^, (17) 

as discussed in [9j [10] it is possible for photons of the opposite chirality (positive if r/ > 0, negative 
if ?7 < 0) if 

|k| > ^- = kth (m^ = 0) (18) 

Here we shall be concerned with the related process e ^ 67 whose kinematics will be discussed 
in detail in the following sections. In fact since we shall be concerned about the possible influence of 
the pseudoscalar background on cosmic ray propagation, which consist mostly of protons, we shall 
discuss the process p — > pj. In this case there is a large hierarchy between the proton momentum 
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and the time variation of the background. We shah see below that the process p — > p7 is perfectly 
possible kinematically. This is a prompt process that is not affected by a slight (in relative terms) 
time variation of the background. The net effect shall not average to zero. 

As for the value of to use, the density of electrons in the universe is expected to be at most 
of the order of Ug = 10"'' cm~'^. Photons pick up a mass of the order of 

Ana—. (19) 

rUe 

This number is of the order of 10~^^ eV, but we shall assume the more conservative limit 10"^^ 
eV, compatible with some astrophysical bounds in the 10~^^ to 10"^"^ eV region [12J. It is a very 
small number, but it is non-zero for sure (it actually depends on the local density of electrons, of 
course) . 



3 Kinematic constraints 

Having found out the different polarizations and dispersion relations in the axion background let 
us now turn to kinematical considerations. Let us consider the process p{p) p{p — k)7(k) with 
= (^E,p),p = IpI and A;^ = (w^^k), k = |k|. Using (|15|) energy conservation leads to 



y/E"^ + /c2 - 2pk cose + Jk"^ + r]k + - E = 0. (20) 



For simplicity we take the positive sign for the polarization, in the understanding that changing 
the sign of r] amounts to exchanging positive and negative chiral polarizations for the photon. 
Let us first consider the case = 0. Then the above energy conservation equation reduces to 

fc2(-4^2 _ pQg2 Q _ _^ k{Ap'^r] + 4771^7?) = 0. (21) 

This equation has the trivial solution k = 0, where nothing happens, and 

4p2 _|_ 4^2 

k = — ^ . (22) 

— 4£'2 -|- 4p2 (,Qg2 ^ _j_ j^2 _|_ 4pj^ 

To find out the kinematical restrictions on k we search for the extrema of the denominator. There is 
only one such extremum which is is cos 9 = —r]/2p. Plugging this into the expression for k one gets 
k = —r]. Since k has to be positive this process is kinematically impossible unless r/ < 0. Of course 
for ?7 = it is also impossible as this is the familiar vacuum case. So the process p{p) — > p(p— k)7(k) 
is actually possible, for positive values of r/ if the emitted photon is of negative chirality. On the 
other hand, if r] is negative, the process becomes possible only for positive chirality photons. Here 
we shall take r] to be always negative, but it is trivial to recover the results for positive r] by 
exchanging chiralities. 



This extremal value (for negative ??) is clearly providing the minimum value for k, kmin = Ivl- 
There is no other extremal point, so the maximum value for k has to lie at one of the two values 
cosO = ±1. It is not difficult to convince oneself (and quite intuitive) that the maximum is found 
for cos 6 = 1 and it corresponds to 

kmax — 13") (23) 
P 

SO kmax ^ E for E > ml/\7]\ and kmax ^ \v\E'^/ml for E < ml/\r]\. 

Let us now consider the case > 0. Conservation of energy leads to the following two roots 
for k 



2E'^ri — 2pm^ cos 9 — rjnii =p 2E^ E-^rf — 4:E'^ml: + + Ap^nv^ cos^ 9 + 2pr]m'^ cos 9 

= . . y [ [ [ [ ^24) 

—AE"^ + cos^ 9 + Apr] cos 9 + rj^ 

The equations above are meaningful only for those values 9 and p providing a positive discriminant 
A. From the condition that A > we get 



where 



P. = ^±^jl-/,. (26) 
—r] —r] y Am^ 

Let us examine these two values p±. Since rj is taken negative, it is clear that p+ > and p- < 0. 
Therefore for sin'^ to be a positive quantity we need 

P>P+=Pth- (27) 

-77 

This is the threshold value below which the process cannot take place kinematically. For p >> pth, 

sin^ 9max ^ ■ (28) 



9max is the value of 9 where the bound in sin 9 of eq. ()25p is saturated. We see that photons are 
emitted in a rather narrow cone. This justifies a posteriori the approximation cos 9 ~ 1 — ^sm'^9 
that has been used. 
At 

9max the square root in ()24p vanishes and — k— — k(^9max) in (|24|) . Then, keeping only 
the leading terms we get 

ki9max) ^—^{1 + —L (29) 

—rj E'^rj —r] 

From eq . (|24|) we work out the value for 9 = 0, which is the minimum value of 9 from the bound in 
eq.(l25]). 



E rj — pmii — E^l E'^rf — Arn^rni^ + 2pr]m'i^ ^2 



T P»Vth 

n — 
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fe+(0) ^ ■ (30) 

+ ^ ^ 2pr]-2ml ^ . 



which is the same result obtained before, and 



E"^-!] — prn^ + E^ E'^rp' — AmpUi'^ + 2pr]m'i^ p»p 

A;_(0) ~ - - 2 ^ (31) 

Ipr] — Irrip —r] 

Now we notice that fc_(0) < k{9max) < ^+(0). To show that kmin = k^{0) and kmax = ^+(0) we 
have to study the derivative of 9 versus k, namely we should have dcos9/dk < for k < k{9max) 
and dcos9/dk > for k > k{9fnax)- 

We isolate cos 9 from the energy conservation relation ()20p 



rnz + rjk — ^E^jrai^ + k"^ + r]k 

cos 9 = — . (32) 

-2pk ^ ' 



and 



d cos 9 _ rn"^ E r]k + 2m'^ p»Pth r/fc + 2m^ 



dk 2k^p 2k^p Lf^ + ^2 + ^2 2k\/7]k + k^ + m^^ 



(33) 



For k k{9 

max) d COS 9 /dk — > and is the only zero which of course means that this value of k 
corresponds to a minimum of cos 9 (i.e. to a maximum of sin 9). On the other hand, for k < k{9max), 
dcos9/dk < and for k > k{9max), dcos9/dk > 0. Then kmin = A;_(0) and kmax = ^+(0) 

We have said before that the actual cone of photon emission is quite narrow for the values of 
and r] relevant in the present discussion, but for the sake of completeness let us discuss the case 
where is small. We use again the expression (|32p and 



2pk'^ -y/ + /c^ + ki] 



The extremum conditions leads to A; = or + A;^ + fcr/ = 0. The first possibility is unphysical 
since it is incompatible with (j32]) . The second one gives 



^- '"^"V^ . (35) 

which is real only as long as 77^ > 4m^. These expressions have a smooth — > limit which has 
already been discussed at the beginning of this section. 



4 Differential decay probability 

Here we shall compute the relevant matrix element and the differential decay probability. The 
calculation is standard but in view of the rather peculiar properties of Chern-Simons QED it is 
probably useful to go into some details. Since the process is possible only for one chirality of the 
photon (depending on the sign of rj) we shall not sum over the final polarization of the radiated 
photon. 
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Using the standard Feynman rules we get 

|M|2 = eh*^e^ti [u{q)u{q)j''u{p)u{p)Y] (36) 

Summing and averaging over the final and initial proton helicities, respectively, and performing the 
trace we get 

\M^ = 2e^-p-k+ [e;{k)e,{k) + e^{k)et{k)] p^^p"} = 2e^ {-p ■ k + 2\e''_{k)p^\^) , (37) 
where e^l(k) is the polarization vector defined in section 2. Recall that [9j 

e^l(k)*e!:(k) + e'l(k)e!:(k)* = -^-^ (38) 

where 

S^"' = {{t] • kf - rfk^)g^"' - (?? • A:)(7?'^F + ri^'k^') + k^r^^'r,'' + rfk^'k". (39) 

Then 

|M|2 = 2e^ {-p-k+ IpP sin^ 9) . (40) 

It is not difficult to see that in the kinematical conditions where the process is possible p ■ k < 0, 
thus guaranteeing positivity. 

The differential decay width will be 

dT{Q) = {27r)^6'^'\q + k-p)^jM^dQ, (41) 

where dQ refers to the final state phase space. The final result is 



1 llfl I Ml^ /T Ikl 1 

dT{Q) = LlL^d\k\ = ^Ll i-p . k + IpP sin2 6)d\k\ (42) 
loTT IpI -C/Wk 2 IpI -C/O^k 

after using the energy-momentum conservation delta and performing the integral over all the angles. 

In the previous expression s'm9 corresponds to 

^"^0= ' -2|p||k| ^ («' 

5 Energy loss and radiation spectrum 

Since in this section no confusion with four-vectors is possible, we revert to the lighter notation 
p = IpI, fc = |k| and w = w^. The relevant quantity to compute is the following 



dE _ 1 

dx V 



dr{Q)w{Q). (44) 



Then, using the previous results, 

j kdk[-^{m'^^ + r]k)+p'^{l-cos'^e)]. (45) 
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dE a 1 
dx 2 p2 



In the previous expression cos 9 is given by ()43p . Then the l.h.s of this identity equals 

- [ ^[-27]k^ - 2mlk^ + Ap^k^ - rj^k^ - mi- AE'^vo^ - 2mlr]k + AEmlu + Ar]Euok]. (46) 



Expanding lo = k\J 1 + ^ + ^ k, k + ^ + recalling the integration limits ^ and ^^^^2 and 
neglecting at this point the effective photon mass we get 

There are two relevant limits 

S < ^ ^ — = f^. 48 

\r]\ ax 4m^ 

At this point it becomes obvious why we have bothered to keep the proton mass all along the 
calculation. 

There are two key scales in this problem. One, discussed in sections 2 and 3, is the threshold 
energy where the process p ^ pj becomes kinematically possible in the presence of a pseudoscalar 
background represented by 7^ 0, namely Eth — 2m^mp/\7]\. The other relevant scale is mp/\r]\ 
where there is a crossover: the energy loss per unit length well below this energy is effectively 
proportional to 7]'^, while well above that scale is proportional to ??. Thus even if we are talking 
about very energetic particles the mass is a relevant parameter when Lorentz violating interactions 
are present. 

If E ^ mp/\r]\ the energy loss is given by 

= exp-^x. (50) 

For T] in the range of values described in section 1, this would give a mean free path in the range 
0(1) to 0(10) kpc. This would imply that cold axions act as a powerful shield againts very energetic 
cosmic rays. Since Eth — 10^^ eV the detection of cosmic rays above that energy would in fact 
impose a rather stringent bound on the combination y/p^/fa- 

However, this is not so because even for the most energetic cosmics, just below the GZK cut-off 
of 10^'^ eV, we are well below the cross-over scale mp/\rj\. In this regime the expression for E{x) is 

E{x) = ^ . (51) 

It is peculiar to see that for extremely large distances E{x) ~ ^, independently of the energy of 
the primary particle. However the distances where this behaviour could be seen are unphysically 



large because is very small. 



Thus the effect of the pseudoscalar background on the reach of cosmic rays is totahy neghgible 
for reasonable energies. Does this mean that the effect is completely invisible? Not quite. 

Let us assume for the sake of definiteness that E{0) ~ 10^° eV. Then, the photon emission 
due to the axion background goes all the way up to Wmax ~ 10~^ eV (for -^/p^/fa just in the 
currently allowed upper limit, and getting smaller as 77 decreases). These energies corresponds to 
radio emission in the submillimetric range (10^^ m or larger). 

For electrons or positrons the numbers are different. The threshold energy for the radiation 
loss due to the interaction with the axion background can be as low as 1 TeV (for ^ 10^^ GeV). 
The radiation emitted in this case would be with wavelengths characteristic of X-ray emission, or 
larger. 

It is not easy to say whether these photons could be detected at all, because even over cosmolog- 
ical distances the rate of radiation is very low for protons (higher for electrons), but the emission 
of cosmic rays is sometimes copious and this could offset the low probability of radiation. The 
shape of the spectrum is completely calculable by the formulae we provide here, and the radiation 
is emitted in a very narrow cone around the charged particle. These two characteristics may be 
instrumental to detect it or, alternatively, to exclude its presence. Excluding it with some confi- 
dence level would place an interesting bound on -^/p^/fa- It should be mentioned that, being low 
energy, this radiation once is produced is not affected by the pseudoscalar background as discussed 
in section 2. 

6 Conclusions and outlook 

In this work we have considered the effect on charged particles of a mildly (compared to the particle 
momentum) time dependent pseudoscalar background. A physical situation worth exploring is the 
influence of a diffuse relic cold axion background on cosmic ray propagation. The effect is model 
independent, universal and can be computed unambiguously. 

The effect is completely calculable in great detail because particle propagation is governed by a 
modification of QED (Maxwell- Chern- Simons Electrodynamics) that is exactly solvable. We have 
determined the kinematical constraints, the characteristics of the emitted radiation and the rate of 
energy loss of charged particles moving in such a medium. Some rather non-intuitive features appear 
and the results, we believe, are interesting per se. We find that for protons the phenomenon may 
appear (depending on the value of ^/fki/fa) at energies around 10'^ TeV, but for those primaries with 
energies that survive the GZK cut-off the suppression brought about by the interaction with this 
'axion shield' is quite negligible, being only relevant for extremely energetic particles. However, the 
'bremmstrahlung' radiation may be measurable. The characteristic energies and extremely strong 
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colimation associated to given primaries may be crucial for that. 

We have also investigated the relevance of this phenomena for electrons and positrons. In this 
case the threshold energy is really around the corner (around 1 TeV, maybe even less as it depends 
on the photon effective mass and on the value of y/p^/ fa) and the radiation could be as energetic 
as X-rays. 

An analogous phenomenon may present itself in a completely different context, namely heavy 
ion collisions at high density and relatively low temperature. As emphasized recently [ISj a phase 
where a parity breaking condensate is present may well appear even at moderate densities. Once 
formed this condensate would be time dependent due to the expansion of the fireball produced 
in the collision. Such a background would have rather non trivial effects on the propagation of 
different particles. It is known to potentially influence the behaviour of scalars and pseudoscalars, 
dramatically changing the spectrum, but it would also influence, via a mechanism similar to the 
one discussed here, photons and leptons. This effect is being investigated [14] . 
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